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Abstract 

Following a similar recent analysis for CP violation in the electroweak sector of the standard 
model, we estimate the naturalness of a magnitude of CP violation (measured by the Jarlskog 
invariant J) close to the observed value in extensions of the standard model with left-right sym- 
metry, such as the Pati-Salam model, where quark mass matrices are not Hermitian in general. 
We construct a simple and natural measure on the space of complex matrices which is both 
geometrically motivated and uses the observed quark mass hierarchy. We find that, unlike in 
the case of the standard model where the observed value for J seemed rather typical, one would 
now expect to observe J < 10~ 7 , clearly in conflict with the observed value J ~ 3 x 10~ 5 . The 
crucial difference in the calculation lies in the non-Hermiticity of mass matrices that modifies 
the measure. We conclude that one would need additional assumptions modifying the measure 
to reproduce the observed value, and that in this sense the standard model is preferred to certain 
classes of left-right symmetric extensions: It does not need additional assumptions to explain 
the magnitude of CP violation. 
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1 Introduction 

In modern theoretical physics one often tries to make statements about "naturalness" or "fine- 
tuning" of the observed values of fundamental parameters, where fine-tuning of a parameter is inter- 
preted as an indication for incompleteness of the theory. Popular examples of fine-tuning problems 
include the quark mass hierarchy and the cosmological constant problem in particle physics. Since 
statements about naturalness are fundamentally of a statistical character, to make them mathemat- 
ically precise one has to assume a well-motivated probability distribution on the parameter space 
relevant for the theory. A fine-tuning problem then indicates that the probability distribution one 
has used should be modified by introducing new physical considerations. As an example, as long as 
observation was consistent with a vanishing cosmological constant, it seemed reasonable to assume 
that a postulated symmetry would constrain it to vanish. With more recent observations indicating 
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that it must be taken to be very small and positive, there seems to be an issue of fine-tuning]]. 
We note that much of the motivation to extend the standard model of particle physics is driven 
by such considerations, that it is by no means necessary to contemplate a "Multiverse" where all 
possible values of a given parameter are actually realised, and that we need not consider anthropic 
arguments. Of course we can only ever observe and make measurements in a single Universe, and if 
a parameter takes a value that appears unlikely maybe this just means that an unlikely possibility 
is realised in our Universe. One is merely doing statistics. But one should recall Bayes' theorem [5] 

_ P(B\A)P(A) 

( 1 ^'E^WW () 

where one can take A as a hypothesis and B as an observation. Once we specify A priori probabilities 
for a full set of possible hypotheses Aj, and compute the probability of observing B given Ai, this 
allows us to make a statement about the probability of the hypothesis following from the observation. 
If the priors P{Ai) are modified, one might get a very different answer for P(A\B). 

The particular example we will investigate in this paper is the apparently weak magnitude of 
CP violation in the electroweak sector of the standard model, measured by the Jarlskog invariant 
J. The observation by Kobayashi and Maskawa [3] that CP violation is only possible for at least 
three quark families has led to a Nobel prize, but the issue of possible fine-tuning in the magnitude 
of CP violation is much less understood. 

It is true that, if one is talking about parameters in quantum field theory, they should really not 
be regarded as constants, but have an evolution with energy scale given by renormalization group 
equations. However, since fine-tuning means a discrepancy of several orders of magnitude, it may 
well be that a fine-tuning problem is present at all energy scales. This is true for the quark mass 
hierarchy [3], and also for the case of CP violation: Recent numerical studies [5] indicate that J 2 
does not run strongly with energy scale, but that the value at extremely high energies (~ 10 15 GeV) 
is merely about twice the value at low energies. It is then meaningful to talk about "naturalness" 
of the value of such a parameter. 

Since any statements one tries to make depend very directly on the choice of measure, it is helpful 
if geometric considerations allow for a natural choice of probability distribution. For instance, if the 
parameter space is a homogeneous space G/H for G a compact Lie group and H a closed subgroup, 
the natural requirement on the measure determining the probability distribution is invariance under 
the left action of G, which leads to a unique measure (up to normalization). Probabilities for a given 
function on G/H to take certain values are then well-defined. 

In the case of CP violation in electroweak theory, the parameter space is the space of Cabibbo- 
Kobayashi-Maskawa (CKM) matrices, a double quotient H\G/H, which makes the problem of 
determining a natural probability distribution more involved than for a homogeneous space G/H. 
This problem was discussed in 6], where several possible choices were investigated. It was found 
that while there is no clearly preferred choice of measure, there always seems to be fine-tuning in 
the observed value for J, unless additional input is used. 



1 For an alternative interesting but presumably non-mainstream viewpoint, see pQ 
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In the second part of [5] (summarised in [7]), the observed values for the quark masses were taken 
into account by considering not the space of CKM matrices, but the space of mass matrices, as the 
fundamental parameter space. This was motivated by the observation that the mass matrices are 
directly linked to the Yukawa couplings and the Higgs vacuum expectation value, whereas the CKM 
matrix is only a derived quantity. The observed values for the quark masses were then taken as 
given, and a probability distribution was constructed that could reproduce these values. The choice 
made for this distribution was as simple as possible in the following sense: The natural group action 
on the space of Hermitian mass matrices is the action of the unitary group U (3) by conjugation, 
AI -> UMW. There is essentially a unique measure invariant under this action. This measure was 
then modified by introducing the simplest possible function that would allow for a modification of 
the expectation values for quark masses fitting observation. No further assumptions were needed. 
It was then found that this simple choice for the measure gave an expectation value for J that was 
remarkably close to the observed value. Hence the conclusion of [7J was that once one assumes the 
quark masses as given, one does not face an additional fine-tuning problem with J. This statement, 
while not new, had been made precise using a geometrically motivated measure on the parameter 
space. 

The calculations done in [5] heavily relied on the fact that the mass matrices in the standard 
model can be taken to be Hermitian without loss of generality. Indeed, the starting point was the 
most natural measure on the space of 3 x 3 Hermitian matrices. However, in left-right symmetric 
extensions of the standard model, such as Pati-Salam, such an assumption can no longer be made 
and the mass matrices have to be regarded as a priori arbitrary complex matricej^]. In this paper, 
we investigate the consequences for statements about naturalness of J by redoing the analysis of 
[6] for general complex matrices. We again use the most symmetric measure on the space of mass 
matrices, here the space of general complex matrices, and modify it in the simplest possible way 
to incorporate the observed quark mass hierarchy. While this is a choice that could of course be 
made very differently, it is a simple choice that uses as few assumptions as possible, and that has 
worked very well for the case of Hermitian mass matrices. What we find is that the resulting 
probability distribution is different, and the result is different too: The observed value for J now 
appears to be unnaturally large, since CP violation should be more heavily suppressed by the quark 
mass hierarchy. One faces a fine-tuning problem, and needs additional assumptions to modify the 
measure appropriately. 

We should point out that the only real input we use from the physical theory (standard model or 
a left-right symmetric extension of it) is, apart from the very definition of CP-violating parameters, 
how the theory restricts the type of mass matrices that appear. In particular, for any extension of 
the standard model (such as left-right symmetric models where parity is the left-right symmetrjlfl) 
that also has Hermitian mass matrices we would not see any modification in the results. 

The structure of the paper is as follows: In section 2, we review CP violation in the electroweak 

2 We thank Ben Allanach for pointing this out. 

3 Thanks to the referee for clarification on this point. 
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sector and detail why mass matrices may be assumed to be Hermitian in the standard model, but 
not if one has an extended left- right symmetry. In section 3, we construct a measure on the space 
of 3 x 3 complex mass matrices, taking into account both the invariance of a natural measure under 
left or right multiplication of a complex matrix by a unitary matrix as well as a decaying function 
that is necessary for the calculation, includes the observed values for the quark masses, and partially 
breaks this invariance. In section 4, we compute the expectation value for the square of the Jarlskog 
invariant J in the given probability distribution, showing that it is much smaller that the observed 
value. This main part is very similar to the calculations in [6] and the reader may benefit from 
comparing with this paper. We summarise in section 5. 



2 CP Violation in the Standard Model and Beyond 

We summarise how CP violation arises, first in the standard model and then in the more general 
case of left-right symmetric extensions, essentially following [5]. In the standard model, the quark 
fields appear as left-handed SU(2) doublets and right-handed SU(2) singlets: 

QjR, q' jR , j = l,2,...,N. (2) 

Here N is the number of quark families, which is arbitrary in the standard model, and normally 
taken to be three. The fields are written in a flavour basis which can be considered unphysical, since 
flavour eigenstates do not correspond to mass eigenstates. 

The coupling of the Higgs doublet H to quarks, through spontaneous symmetry breaking by the 
Higgs potential, gives masses to the quark fields: 

N 

£m gS s ^ - ( m jkqjLlkR + m' jk q' jL q' kR J + h.c. (3) 

The mass matrices m and to' are determined by the original Yukawa couplings and the Higgs vacuum 
expectation value. Thus, it seems appropriate to regard either the set of Yukawa couplings together 
with the Higgs vacuum expectation value or the collection of elements of to and to' as fundamental 
parameters of the theory. This viewpoint was supported by the results of [6] . 

This part of the Lagrangian is formally CP invariant if and only if to and to', which are so far 
arbitrary complex matrices, are real. Since this condition is not satisfied in Nature, one has formal 
CP violation. However, as remarked before, the Lagrangian has been written in the unphysical 
flavour basis. One can, for general to and to', pass to a different basis, namely the basis of mass 
eigenstates, by diagonalising the mass matrices with unitary matrices: 

m=UlAU R , m' = (U' L )^'U R , (4) 

thus the basis of mass eigenstates is related to the previously considered basis by 

<7Lphy S = U L q L , tephys = U R q R , etc. (5) 
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It is always possible to choose the unitary matrices so that A and A' are real, and in this new 
basis this part of the Lagrangian is invariant under C and P separately, and hence also under CP. 
However, the electroweak Lagrangian also contains charged current terms mixing up- and down-type 
quarks, coupled to the W boson fields via (we are now using the basis of mass eigenstates) 



where V := Ul{U' l )^ is the Cabibbo-Kobayashi-Maskawa (CKM) matrix. In the basis of mass 
eigenstates, this term Xq is not invariant under CP unless V is real. Since we consider the mass 
eigenstates as physical, CP is violated through these charged current terms. 

An important observation made in [5] is that one can redefine the right-handed quark bases by 
arbitrary unitary transformations, 



obtaining a new basis which is to be regarded as equally physical. This is due to the absence of 
charged current terms involving right-handed quarks, since they are singlets under SU(2). It is 
therefore no loss of generality to set Ur — Ul and U' R — U' L in (Q}, and to assume that m and m! 
are Hermitian. 

A natural way to extend the standard model is to assume the existence of a second SU (2) 
symmetry which acts on the right-handed quarks, as in the Pati-Salam model jTO] . In such extensions, 
one adds a term © for right-handed quarks to the Lagrangian. This has the important consequence 
that a general transformation ([7J for arbitrary unitary transformations O and O' can no longer be 
regarded as giving an equivalent quark basis, since it modifies this new charged current term. The 
mass matrices cannot in general be taken to be Hermitian, but are arbitrary complex matrices. 

There are now also two possibly CP violating terms, and two CKM matrices. We focus on 
V = Ul(U' l Y which involves the processes that are actually observed and disregard Vr = Ur{U' r )^ 
in the following. 

Mathematically, V is an element of SU(N), but since the phases of the quark fields are arbitrary 
(even for non-Hermitian mass matrices), V is only defined up to left or right multiplication by a 
diagonal element of SU(N), i.e. an element of the maximal torus [/(l)^ -1 . The space of CKM 
matrices is therefore the double quotient U(l) N ~ 1 \SU(N)/U(l) N ~ 1 , characterised by (N — l) 2 
parameters, out of which ^N(N — 1) may be taken to be real (Euler) angles and the remaining 
i(iV — l)(N — 2) appear as complex phases. It follows that the matrix V can be taken to be real 
for N = 2, and so in the formalism explained here one needs at least three quark families to have a 
possibility of CP violation. Kobayashi and Maskawa were awarded the 2008 Nobel Prize for using 
this observation to predict the existence of a third quark family (3J. We set N = 3 in what follows. 

The mathematical theory of observable measures of CP violation in the standard model was 
developed by Jarlskog [8]. She showed [TT] that all necessary and sufficient conditions for CP 



fd L \ 



Xc ■■= (Wl - iWt)J£ + h.c. , J? := (u, c, t) L -fV s L 

w 



(6) 



U R -> OU R , U' R -> 0'U' h 



R ■ 



(7) 
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violation can be summarised as the following condition on the commutator of the Hcrmitian matrices 
m, m': 

detC:=det(-i[rn,m']) ^0. (8) 

One finds that explicitly 

det C = -2J(m t - m c )(m c - m„)(m„ - m t )(m b - m s )(m s - m d )(m d - m b ) , (9) 

where J :— 3tn(Vii ^22^1*2^2*1) is the Jarlskog invariant which is invariant under left or right multi- 
plication of V by a diagonal matrix, i.e. an element of U(l) 2 . The geometrical interpretation of the 
quantity J is given by the so-called unitarity triangles. These express the requirement on V to be 
unitary, so that for instance 

{W*) 12 = v n v 2 \ + V 12 V* 2 + V 13 V* 3 = . (10) 

In the complex plane, the three complex numbers that sum to zero form the sides of a triangle. The 
absolute value \J\ is twice the area of this triangle. Since there are different unitarity triangles cor- 
responding to different elements of V V ^ , all with the same area, there are several ways of expressing 
J in terms of the elements of V. A general formula is given by 

J e *M e jki = ?m{V aj Vp k V* k V^) . (11) 

The quantities describing the CKM matrix which are invariant under rephasing of the quark fields 
are J and the absolute values \V a j\. 

In the general case where m and m' are not assumed to be Hermitian, the corresponding quantity 
would be 

det C := det (— i [mm\ m! W^] ) = -2J(m(-ra^)(mj-m^)(m^-m()(mj-mj)(mj-m^)(mj-mj) , 

(12) 

which of course leads to the same conditions on the mass matrices as before (given that we took all 
masses to be positive before). In the literature, the use of C is perhaps more common than the use 
of C, and one may well argue that this second measure of CP violation should be considered more 
fundamental as its value does not depend on the arbitrary signs of the mass terms in the Lagrangian. 
One common parametrization of the CKM matrix, which we use in the following, is given by 

^ cos y cos z cosy sin z e~'™siny\ 

V= — cos x sin z — e lw sin x sin y cos z cos x cos z — e lw sin x sin y sin z sin x cosy , (13) 
\ sin x sin z — e™ cos x sin y cos z — sin x cos z — e™ cos x sin y sin z cos x cos y J 

where the ranges of the Euler angles x,y,z and the complex phase w are 

0<x,y,z<^, 0<w<2ir. (14) 
An arbitrary SU(3) matrix can then be written as 

U = T L VT R , (15) 



G 



where 

T L = diag(e 2i?) , e^^, e -i ^) , T fl = diag(e i(r+t) , e i(r ~ 4) , e~ 2ir ) , (16) 

and p, q, r, and t are phases which can take all values between and 2tt. This shows that the coordi- 
nates x,y,z, and w indeed parametrise representatives of the double quotient U(1) 2 \SU(3)/U(1) 2 . 
In this parametrization, the Jarlskog invariant J is given by 

J = — sin 2.t sin 2z sin y cos 2 y sin w . (17) 

It appears that the observed value for J is very small, since the maximal value would be « 0.1, 
whereas in Nature [H] 

J = 3.05±g;|g x 10~ 5 . (18) 

In a general discussion where the values of the quark masses are not fixed, J is not an appropriate 
measure of CP violation, since even with nonvanishing J one could have CP conservation if, for 
example, the up and charm quark masses were coinciding. It was suggested in |13] to use an 
appropriately normalised form of det C, namely 

a CP = 3V6 (T ^ 2 ^ 3/2 (19) 

for three quark families as the unique basis independent measure of CP violation. This is a dimen- 
sionless number which takes values between —1 and +1, and is again observed to be very close to 
zero. When written out in terms of the CKM matrix parameters and quark masses, it is a rather 
complicated expression which is therefore not extremely useful in practical computations. In the 
present analysis, we assume the quark masses as known and regard J as the measure of CP violation. 



3 Measure on the Space of 3 x 3 Complex Matrices 

In this section we determine a natural measure on the space of 3 x 3 complex matrices in order 
to make statements about likely or natural values for the magnitude of the Jarlskog invariant J. 
Following 6 , this measure is a product of the geometrically most natural measure with maximal 
symmetry and a factor that involves the observed values for the quark masses. This second factor is 
introduced for two reasons; firstly, to make the total volume of the parameter space finite, secondly, 
to allow for a modification of the distribution on the space of CKM matrices through the quark 
masses. The analysis is analogous to the case of 3 x 3 Hermitian mass matrices discussed in [BJ, and 
we will see shortly where differences arise that eventually lead to different results. 

The space of 3 x 3 complex matrices has a natural left and right action by U (3) corresponding 
to changes of basis. We first determine a metric invariant under these group actions. We start with 
Jarlskog's representation of an arbitrary complex matrix as 

M = UlDU R , (20) 

where Ul, Ur are unitary and D = diag(Di, D2, -D3) is real diagonal. Here, as suggested in [TT], M 
is the dimensionless mass matrix M = m/A, where A is a scale which may be chosen for convenience. 
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A natural choice would be A = m t for the up-type or A' = mj for the down-type quarks, but since 
we in principle allow arbitrary values for the quark masses we leave A arbitrary. Clearly, Ul and 
Ur are only defined up to simultaneous left multiplication by a diagonal unitary matrix 

U L ^AU L , Ur^AUr, AeU(l) 3 (21) 

which reduces the number of (real) parameters from 21 to 18. There are additional discrete ambi- 
guities, corresponding to a permutation or change of sign of the elements of D, given by elements of 
the group S3 x Z 3 ,, where 63 is the symmetric group of three elements. (We will see shortly that the 
measure vanishes whenever elements of D coincide up to sign, so we can restrict to matrices with 
D\ ^ D\ 7^ D\ 7^ D\.) Hence we can identify the relevant subspace of the space of 3 x 3 complex 
matrices with R 3 x {U{lf x 63 x Z|)\(17(3) x £7(3)) ~ R 3 x (U(l) 3 x 6 3 )\(t/(3) x 17(3)). However, 
since all expressions will only involve absolute values of the elements of D and D' (precisely due to 
the Z 3 , symmetry) , we will integrate over all of R 3 for simplicity. This only leads to an extra factor 
8 which drops out in expectation values. The metric we will use to determine a measure is 

ds 2 = Tr(dM • dAf+) , (22) 

which clearly is invariant under M — > OMO 1 for O, O' e U(3). To evaluate this, use 

dM = -UldU L ulDU R + uldDU R + ulDdU R , 

dikft = -U R dU R U R DU L + U R dDU L + U R DdU L (23) 
and [D, dD] = to obtain 

Tr(dM-dMi) = Tr(dD-dD)-Tr(D 2 (dU L Ul) 2 )-Tr(D 2 (dU R U R ) 2 ) + 2Ti(DdU L UlDdU R U R ) . (24) 
We introduce right-invariant one-forms 

dU L Ul=i\ a T a L , dU R U R = i\ b T b R , (25) 
where Ai, . . . , As are the Gell-Mann matrices, and 



A 9 = 



(1 





o\ 





1 





V° 





V 



(26) 



so that iA a are a basis for the Lie algebra u(3). Then 

ds 2 = Ti(dM-dM^) - Tr(d7J.d7J) + Tr(7J 2 A (a A b) )(T^T^ + r^)-Tr(7JA a 7JA b )(T£4 + T^) . (27) 
The only nonvanishing traces are 

Tr^AiAi) = Tr(7J 2 A 2 A 2 ) = Tr(7J 2 A 3 A 3 ) = Tr(7JA3^A 3 ) = D\ + D\ , 
Tr(7J 2 A (3 A 8) ) = Tr(7JA (3 7JA 8) ) - ±(D 2 -D 2 ) , Tr(7J 2 A (3 A 9) ) = Tr(7JA (3 7JA 9) ) = \[\{D\-D 



I), 
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A2 

Tr(£ 2 A (8 A 9) ) = TV(£>A (8J DA 9) ) = ^-(D 2 + D 2 - 2D 2 3 ) , 
Tr(D 2 \ 4 \ 4 ) = Tr(D 2 \ 5 \ 5 ) = D\ + D 2 , Tr(L> 2 A 6 A 6 ) = Tr(D 2 A 7 A 7 ) = D\ + D 2 3 , 
Tr(Z? 2 A 8 A 8 ) = Tr(DX 8 DX s ) = ^(D 2 +D 2 +4D 2 ) , Tr(£ 2 A 9 A 9 ) = Tr(D\ 9 D\ 9 ) = ^(D 2 +D 2 +D 2 ) , 

Tt(DXiDXi) = Tr(DX 2 D\ 2 ) = 2D Y D 2 , 
Tr(£ 2 A 4 A 4 ) = Tr(£ 2 A 5 A 5 ) = 2D 1 D 3 , Tr(£ 2 A 6 A 6 ) = Tr(Z? 2 A 7 A 7 ) = 2D 2 D 3 . (28) 
The metric can be written in the form 

ds 2 = dDl+dDl+dDl + ^D.-D^irl+T^f + ^ + D^irl-r],) 2 

+\ [( Dl - D 2 ) 2 (r 2 + r 2 R ) 2 + (A + D 2 ) 2 (t 2 - r 2 ) 2 
+( Dl - D 3 ) 2 (4 + 4) 2 + ( Dl + D 3 ) 2 (ri - r^) 2 ] 
+~ [{Dx - D 3 ) 2 {rt + rl) 2 + (D 1 + D 3 ) 2 (r 5 L - ) 2 
+(D 2 D 3 ) 2 (4 + 4) 2 + (D 2 + D 3 ) 2 {rl - r%) 2 } 
+ \{D 2 D 3 ) 2 {rl + rlf + \{D 2 + D 3 ) 2 (rl - r?) 2 

- 4) + ±(4 - 4) + -L(ri - 4) 

r l) + -L(r!-rI)J . (29) 

The fact that the metric only depends on r 3 — rj^ etc., and not on r 3 + rjj etc., again reflects the 
f/(l) 3 that has to be factored out. As is easy to show, the volume form is proportional to 

(D 2 - D 2 ) 2 (D 2 - D 2 ) 2 (D 2 - Dl) 2 \DxD 2 D 3 \ dD 1 AdD 2 AdD 3 At* A A . . . A (rg - r|) A (r| - t%) . 

Note that this expression only depends on the absolute values of Di, as expected. Since the range of 
the Di is infinite, integration over these coordinates will give an infinity, so that a function decaying 
sufficiently fast for large \Di\ is introduced. A natural and simple choice is a Gaussian. 

Since there are actually two integrations over the space of mass matrices, corresponding to 
two mass matrices for the up-type and down-type quarks, the general expression for the measure 
considered in [6l [7] , where M was assumed to be Hermitian, was 

DM DM' exp(-Tr(M 2 A))exp(-Tr((Af') 2 ^)) , 

where DM was the natural measure on the space of 3 x 3 Hermitian matrices, with A and A' 
Hermitian with non-negative eigenvalues and commuting. A and A' could then be diagonalized 
by the transformation U — > UW, U' — > U'W without changing any measurable quantities. The 
Gaussian broke the symmetry of the measure from invariance under (U(3) x U(3)) to invariance 
under ([/(f) 3 x U(l) 3 ). 
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The analogous procedure in the case of general complex mass matrices is to use the measure 

DM DM' exp(-Tr(A/A/U + M'(M' ) t A')) = DM DM' eiq>(-Tr(ulD 2 U L A + {U' L )\D') 2 U' L A')) , 

so that the measure is still invariant under the right C/(3) action on complex matrices, but the 
invariance under the left f/(3) action is again broken to U(l) 3 . This symmetry breaking is necessary 
to obtain a distribution that can reproduce the observed quark masses. As before, we assume A 
and A' to be Hermitian with non-negative eigenvalues and commuting and use the redefinitions 
U L -> U L W L , U' L -> U' L W L to diagonalise A and A'. 

Since the important measurable quantity involving the unitary matrices is the CKM matrix 
V = UlU'l, everything is independent of the Ur parameters and we integrate over these, obtaining 
a constant which is irrelevant in the averaging process. 

We are left with integrating over the space of possible Ul, the coset t/(l) 3 \£/(3) = U(1) 2 \SU(3) 
(both the Gaussian and the CKM matrix are invariant under left multiplication of Ul by an element 
of C/(l) 3 ), and the volume form is proportional to 

(Df - D\f{D\ - Dlf{D 2 2 - Dfj^UiDaDa] dD 1 A dD 2 A dD a A r\ A t 2 A t a l A rf A r\ A r\ . 

In terms of the coordinates on U(1) 2 \SU(3) introduced in section 2, the wedge product of right- 
invariant forms gives the usual bi- invariant measure on SU(3), so that we finally get 

DM = ( Y[i D i ~ D jf 1 \DiD 2 D 3 \ sin2x cos 3 y amy sin2z dD x dD 2 dD 3 dx dy dz dw dr dt. (30) 

We also have to take the discrete ©3 symmetry into account, integrating only over one sixth of 
the homogeneous space U(1) 2 \SU(3), corresponding to 

7T 

< y < arctan(sin x) , < x < — . (31) 

This restriction amounts to removing unitary matrices that permute the elements of D and hence 
to fixing an ordering. 

Comparing to the case of Hermitian matrices [7] , the measure involves higher powers of elements 
of D. Therefore, we would expect a stronger influence of the quark mass hierarchy on expectation 
values of (powers of) J, namely a stronger suppression of large values of J. This naive expectation 
will be confirmed in the next section. 



4 Results for J 

The calculations go through exactly as in the case of the measure used in [6] , and our analysis will be 
completely analogous. We will obtain analytical approximations of the relevant integrations whose 
validity is then verified by numerical integration (using Mathematica). 

We want to compute the expectation value for J 2 in the given probability distribution (all odd 
powers of J average to zero), 

_ J DM DM' e~^ MMtA '>- Tr( - M ' ( - M "> fA ">J 2 (M,M') 

\ J / - J DMDM , e -Tr(MMtA)-Tr(M>(M')iA>) ' ( 32 ) 
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where DM is given in (|30p and DM' is the same expression in terms of primed variables. It should 
be clear that multiplying A or A' by a (non-negative) constant is the same as rescaling the arbitrary 
scales A and A' used in defining Di or D[. Hence we may assume, without loss of generality, the 
form 

(\ \ (l \ 



.4 







A' 



/i" 2 
\0 







(33) 



-2 



Mo 2 
\0 

We have introduced four dimensionless parameters /z c ,/z u ,/i s , and [id that we are free to choose so 
as to reproduce the observed values of the quark masses as expectation values in our distribution!*) 
On dimensional grounds alone, we expect that we will have to set fi c ~ m c /m t etc. In particular, 
we will have [i u <C 1 and fid -C 1. The precise relation between the fi parameters and expectation 
values for (squared) quark masses will be determined shortly. 

Now J 2 , written in terms of coordinates on M and M' , is a very complicated expression. However, 
since 



Tr(MMU) = £ DfAAULll = sin 2 y + .. 

ij 



Tr(M'(M')^A') = ^f- 



M 2 d 



sin 2 y' 



(34) 



where we assume (i„«l and ^ < 1 and Aj denote the diagonal elements of A, the integrand in 
the numerator and denominator of (j3"2")l is negligibly small unless y and y' are close to zero. We 
therefore use the approximation smy rj y in the measure DM, setting y = y' = in the remaining 
part of the integrand. The integration over y and y' can be easily performed, and everything is 
independent of w and w' which may hence be dropped as integration variables. One is left with the 
integral 

J dD dD' J d 4 x dV sin 2x sin 2z sin 2x' sin 2z' (^-^i^UAU^ +(d') 2 u' a'u'*) jj/^ 



(J' 



where 



\y=y'=0 



fdDdD' f d 4 xd 4 x' sm2xsm2zsm2x' sm2z' (e- T *( D2UAUt +( D ') 2u ' A ' u '^) I , „ 



tt/4 tt/2 2tt 2tt 

d A x = dx dz dr dt 



oooo 



(35) 
(36) 

(37) 



and similarly for J d 4 x' . The shorthand dD denotes the measure over the Di, namely 

dD= ^J[{Df-D*) 2 ^j \DiDvDz\dD x dD 2 dDs. 
It is possible to analytically integrate over both copies of K 3 in (|55"j) , using 

oo oo oo 

W := / dD x J dD 2 J dD, {D\ - Dl)\D\ - Dl)\Dl - D\f\D x D 2 D,\e-^-^-^l 



24 

757575 (2(Ci 2 ^2 2 + e^i + zizMi + e 2 + e 3 - ^ - ^ - ^ 
-aeik&fe 3 + e 2 + $ - - tUi - - - - e 2 2 a) - 



(38) 



4 We emphasise that we do not make any predictions about quark masses, but use them as an input to modify the 
probability distribution on the space of all mass matrices. 



11 



The explicit expression for J at y = y' = is 



1 



J ( U > U %= y >=o = 2*2**2*' {4,8%, sin(3f + t) + clc z ,s 2 z , sin(3f - t) 



-c z s z s z >(c z , [sin(3f + 1) + sin(3r — 3i)] — s 2 z , sin(3f + i)) 



+c z c z ,s 2 z (c 2 z , sin(3f - t) - s% [sin(3f + 3f) + sin(3f - (39) 

where s x = sin a;, c z > = cosz', etc., f = r — r', and t = t — t' . Integrating (l39l) over r,r',t, and t' 
indeed gives zero, which is why we choose to use J 2 . 

To fix the parameters appearing in the matrices A and A', we first observe that expectation 
values for squared mass matrices, again in the approximation y 1, take the relatively simple form 

7T/4 7T-/2 



L 3 dD D\ J dx J dz sin 2x sin 2z ( e 



,-Tr(Z> 2 (/yl(7 t ) 



l v =0 



(Df) « ^ ^ . (40) 

N 11 tt/4 tt/2 

o o 



ly=0 



The denominator is explicitly 

tt/4 tt/2 



Id '■= J dx J dz sin 2x sin 2z /c^^ , (41) 



where 



£i = Ai cos 2 z + A2 sin 2 z , £2 = Ai cos 2 x sin 2 z + A2 cos 2 x cos 2 z + A3 sin 2 x , 

£3 = A\ sin 2 x sin 2 z + A 2 sin 2 x cos 2 z + A3 cos 2 x , (42) 

with A3 3> A 2 3> Ai. The integral is dominated by very small x and z (we cannot have x = -|), 
and we can approximate Id well by only keeping the terms of leading order in x and z in the 
trigonometric functions, and 

Axfofe * ^5 >< « (Al +^z2 ) 5 ( 4 l 2+ A 3 X 2 )3A3 ' (43) 

which are the leading terms (as we shall see, the combination A3X 2 is effectively of order A2 etc.): 

7T/4 7T/2 

/jD ~ ~I~ J dx J dzAxz ^ Al + A 2 z 2 y 5 (A 2 + A3X 2 )- 3 


00 00 

« ^-JdxfdZ(A 1 + A 2 Z)- 5 (A 2 + A3X)- 3 




48 1 1 6 



A 3 AA 2 A\ 2A 3 A 2 A\A\A\ 



(44) 



This is very well reproduced by numerical calculations. Similarly, we find 

00 00 

Id(D 2 ) « ^JdX J dZ(A 1+ A 2 Z)- e (A 2 +A3XY 





240 1 1 24 



A 3 5A 2 Af 2A 3 A 2 A? A^A 2 ' 
12 



(45) 



hence 



<^)«-±. (46) 



Redoing the same calculation for D2 and D3 gives 

Using the form (j33|) for A, this is 

(U?)«4, (# 2 2 )«2 M 2 , (Dl)n»l. (48) 

If (Df ) is to reproduce the squared top quark mass in units of A, our reference scale for the up-type 
quarks must be A = hint- Then setting m 2 /A 2 = 2/i 2 determines fi c , etc. Apart from numerical 
prefactors of order one, the fi parameters indeed correspond to the quark masses one wants to 
reproduce in the probability distribution. Note that factoring out the ©3 above corresponds to 
fixing an ordering of the quark masses, so that it is possible to obtain unequal expectation values for 
D\ , D2 and D3. Integrating over the whole of U (1) 2 \SU(3) would mean to also include permutations, 
so that necessarily (Df) = (£)|) = (D 2 ). 

Again, because we have to consider the dependence of masses on the energy scale in quantum 
field theory, described by the renormalization group, there is some ambiguity in what is meant by 
the "quark masses" we want to reproduce. Following [14] . for example, we take all the quark masses 
evolved to the scale of the Z boson mass. These are given in [4]: 

(m„, m c , m t ) = (1.27±g;|§ MeV, 0.619 ± 0.084 GeV, 171.7 ± 3.0 GeV) ; 

{m d , m s ,m b ) = (2.90ti; 2 g MeV ' 55 ±i5 MeV > 2.89 ± 0.09 GeV) . (49) 

We use the central values 

{m u ,m c ,m t ) := (1.27 MeV, 0.619 GeV, 171.7 GeV) ; 

(m d , m s ,m b ) := (2.9 MeV, 55 MeV, 2.89 GeV) . (50) 

The mass scales A and A' are now fixed by setting {D\) — (m t /A) 2 and ((Z^) 2 ) = (m& /A') 2 . By 
comparing the results obtained by numerical integration with the values we want to reproduce, we 
can then fix the parameters /i c , fi u , fj, s and /id- 

In the case of the positively charged top, charm and up quarks, which exhibit a more extreme 
quark mass hierarchy, we find that numerical calculations (using Mathematica) reproduce the 
results we have obtained analytically very well. For the negatively charged quarks, we find numer- 
ically that we have to use different relative factors to reproduce the observed masses. Comparing 
the numerical results with (|50[) . we fix the parameters appearing in A and A' to 

5 i 2 (—\ «2.60xl0~ 5 , ul=4(— y ] «2.19xl0~ 10 , 
\m t J \m t J 

£ = (rn£\ Rj362xl0 -4 ) 2 4 (HU\ a 4.03X10-*. (51) 
\m b J \m b J 
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In order to obtain an analytical expression for expectation value of J 2 , the next approximation is 
that the main contribution to the integral (I35[) will come from small z. This again is seen by writing 
out Tr(MAfU) and using the mass hierarchy. We only take the term in (|39p that is non-zero at 
z = into account, setting sin2z « 2z in the measure. 

Averaging over r, i,r', and t' gives a factor of 1/2, as one might have expected, and therefore we 

use 

J small z '■— 2 sm xcos xsvo. x cos x cos z sin z (52) 

for our calculations. Within this approximation for J, still taking f^^ 2 ^ 3 ~ 48^j~ 5 ^ 3 ^ 1 , the 
numerator of (|35p is the product (using again that only small z contributes) 



tt/2 tt/2 

2z f , i sin 2z' cos 2 z' sin 4 z' 



1152 X J dz— — ; 777"777 x / ,/.' 



rr/4 

x / dx 



(A x + A 2 z 2 f/ 2 J Mi cos 2 z' + A 2 sin 2 z') 


sin 2x sin 2 x cos 2 a; 





tt/4 



(A 2 cos 2 x + ^3 sin x) 3 (A3 cos 2 x + A 2 sin x) 



sin 2x' sin 2 x' cos 2 x' 



x / dx 5 5 (53) 

(Aj cos 2 x' + A' 3 sin x') 3 (^ 3 cos 2 x' + A' 2 sin x') 

The first two factors are l/(4AfA 2 ) and 1/(12(A' 1 ) 2 (A' 2 ) 3 ), respectively; for the other two (which 
have the same form) we change variables to X — cos 2 x to obtain 



(A 2 X + A 3 (l - X))3(A 3 X + A 2 (1- X)) 2A\A 2 



(l-^ + 0((A 2 /A 3 ) 2 )) ; (54) 



1/2 

the expressions for the denominator are similar but simpler. Putting everything together, we obtain 
the approximation to lowest order in quark mass ratios 

/T 2 v „ l ^K) 2 _ 4 m» 2 

Wsmall z/ ~ g " 3 m 4 m 2 » 

where the numerical factors appearing in the last line come from the different factors chosen in (|51l) . 
Note that the top quark mass does not appear in this approximate result. This compares with the 
scaling behaviour obtained in [5] (for Hcrmitian mass matrices), 

2 \ m 2 m u m d 

(-'small z) ~ m 3 mc • W 

For numerical calculations we use both the simplified expression J 2 ma u 2 and the expression for 
J given in (|39[) . We find that for the first quantity, the numerically evaluated expectation value, 
(-'small 2) ~ !- 89 x 10 -15 , is about 94% of (|55]l. and the numerical result for (J 2 ) (taken at y = t/ = 0) 
is 

(J 2 ) w 2.07 x 1CP 15 , (57) 

which gives 



A J = \J(J 2 ) ~ 4.55 x 1CT 8 (58) 
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which is now almost three orders of magnitude smaller than the observed value (|18p . Assuming a 
Gaussian distribution peaked at zero, we get 

P(\J\ < 1(T 7 )«97%, (59) 

When the measure presented here is used, there seems to be extreme fine-tuning in J, but now we 
would say that one observes unnaturally large CP violation! This result may look surprising, given 
that the maximal value for J is around 0.1 and the observed value just 3 x 10 -5 , but it shows how 
strongly the quark mass hierarchy suppresses large values of J in our distribution. 

5 Summary and Discussion 

We have tried to estimate the naturalness of the observed magnitude of CP violation in the elec- 
troweak theory under the assumption that there is a left-right symmetry which implies that the 
quark mass matrices can not in general be taken to be Hermitian. We have constructed a proba- 
bility distribution on the space of 3 x 3 complex matrices which takes into account the geometrical 
structure of this space, but also includes a Gaussian factor which makes the total volume finite 
and leads to expectation values for quark masses that can reproduce the observed values if four 
free parameters are fitted accordingly. While this is a choice we make, and all results depend on 
this choice, our measure is the combination of a maximally symmetric measure, invariant under a 
redefinition of a complex matrix by left or right multiplication by unitary matrices, and a Gaussian 
incorporating the observed values of the quark masses. We would have to make additional rather 
strong assumptions to motivate a different choice of measure that would differ appreciably from this 
simple construction. It may well be that such assumptions are justified by the underlying mecha- 
nism determining the mass matrices, but we do not know of such a mechanism yet. Furthermore, 
such a simple choice for the measure was shown in [6], where only Hermitian mass matrices were 
considered, to lead to expectation values for the Jarlskog invariant J that make the observed value 
appear typical. 

The conclusion for general complex mass matrices, as shown here, is rather different. Using the 
given probability distribution, one would now expect J to be about three orders of magnitude smaller 
than the observed value. Hence, there is a fine-tuning problem: Without further assumptions, 
a fundamental theory leading to a left-right symmetric electroweak sector at low energy should 
generically be expected to reproduce very weak CP violation. Invoking the principle of Occam's 
razor, "entia non sunt multiplicanda praeter necessitatem," we would like to conclude that, 
only looking at possible explanations for the magnitude of CP violation in the electroweak sector, 
the standard model should be preferred to left-right symmetric extensions such as Pati-Salam: In 
the latter one needs additional assumptions on the fundamental parameters that resolve the issue of 
observing "unnaturally large" CP violation, that are not necessary in the standard model, or any 
extension of it that allows a restriction to Hermitian mass matrices only. 

Although we have tried to argue that our results are independent of renormalization group flow 
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since the relevant quantities do not run strongly with energy scale, there is another subtle issue: 
The low-energy limit of a left-right symmetric extension with non-Hermitian mass matrices would 
still be the standard model, where mass matrices can be assumed to be Hermitian, leading us back 
to the measure considered in [6l [7]- It would be desirable to incorporate this dependency of the 
assumptions one has used to construct the measure on energy scale into the analysis, namely to 
use a measure which depends on energy scale also. A starting point would be a quantification of 
"non-Hermiticity" that could then flow from zero at low energies to some non-zero value at high 
energies. At present these ideas are however somewhat vague, so that we will have to leave them to 
exploration in future work. 
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